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The  modern  theory  of  diffusion  and  precipitation  of  a  heavy  ingre¬ 
dient  considers  only  gravitational  precipitation.  However,  the  wind, 
whose  speed  u  can  be  tens  and  hundreds  of  times  greater  than  the  velo¬ 
city  w  of  the  gravitational  precipitation  of  the  particles,  blows  both 
within  and  above  the  covering  of  vegetation.  Hence,  both  gravitational 
an'd  inertial  precipitation  must  play  an  important  role  as  far  as  coarsely 
dispersed  aerosols  are  concerned. 

In  field  tests  of  precipitation  of  coarsely  dispersed  aerosols1  to 
determine  the  role  of  inertial  precipitation,  20  wind-vane  collectors 
were  set  up  on  the  test  area.  Each  of  them  was  fitted  with  a  horizontal 
glass  plate  and  a  vertical  slide.  The  drops  were  intended  to  fall  on 
the  surface  of  the  horizontal  plate  by  gravitation,  and  it  was  expected 
that  any  precipitation  on  the  windward  side  of  the  vertical  glass  slide 
would  be  the  result  of  inertia.  The  slides  of  the  wind  vanes  were 
mounted  near  the  upper  end  of  a  relatively  sparse  growth  of  vegetation, 
with  the  average  height  of  the  plants  being  n  =  30  cm.  On  the  basis  of 
the  results  of  a  microscopic  examination  of  the  slides  for  each  aerosol 
fraction,  it  was  possible  to  determine  the  mean  value  of  the  ratio 
nv/nyl  for  20  points,  where  ny  is  the  average  number  of  drops  of  a  given 

fraction  which  land  on  a  unit  area  of  the  windward  side  of  the  vertical 
slide,  while  is  the  same  value  for  the  top  side  of  the  horizontal 

slide.  The  values  of  n  ,/n,  ,  obtained  for  fractions  with  different 

v  hr 

average  drop  diameters  d  at  various  wind  speeds  u(h)  at  a  height  of  30  cm, 
are  listed  in  Table  1. 

-  It  is  clear  from  the  table  that  the  density  of  the  deposits  of 
drops  of  the  same  type  was  several  times  greater  on  the  vertical  slides 
(which  can  be  viewed  as  crude  models  of  leaves)  than  on  the  horizontal 
ones  >  !)•  This  result  is  a  direct  indication  that  inertial  pre¬ 

cipitation  of  an  aerosol  on  vegetacion  can  prevail  over  gravitational 
precipitation. 

T'Conducted  with  the  participation  of  I.  F.  Yevdokimov,  V.  M.  Krasil'nikov, 
K.  P.  Mikulin  and  Z.  M.  Yuzhniy. 


Table  1 


d.  V- 

Mtrs. 
m/sec  > 

13  . 

25 

42 

58 

75 

52 

108 

133 

107 

217 

203 

1,8 

0,84 

2,8 

4,8 

4,4 

4,5 

4,3 

4,1 

•>,6 

2,5 

2,1 

2,2 

3,6 

4,6 

8,3 

7,1 

7,5 

5,0 

3,2 

2,4 

3,0 

3,1 

2,9 

4,1 

1.6 

•2,6 

5,2 

6,1 

5,6 

0}0 

4,4 

3,2 

2,4 

Oy  O 

1.2 

0,76 

The  equation  for  the  material  balance  of  the  additive  for  an  ele¬ 
ment  of  the  plant  layer  having  a  length  dx  and  a  height  h  is 

| C(x,z)un(z)dS  z)dS  +  \P{x,  z)dV  =  0, 

where  C  is  the  concentration  of  the  additive;  N  is  the  diffusion  current 
of  the  additive;  P  =  - [a(z)&(z)u(z)  +  z)w\  C  (x3z)  is  the  loss  of 

the  additive  within  the  plant  layer  due  to  inertial  and  gravitational 
precipitation;  a  is  the  coefficient  of  entrapment  of  the  particles  by 
the  plants;  &  is  the  specific  area  of  the  projection  of  the  plants  on  a 
plane  normal  to  u;  0^  is  the  specific  area  of  the  horizontal  projection 

of  the  plants. 

Following  integration  with  averaging  over  z  and  several  simplifica¬ 
tions,  we  obtain  the  following  regional  condition  at  the  upper  limit  of 
the  cover  of  vegetation  z  -  h: 


dc  ( x ,  h)!dz  —  nC  [x,  h), 

where 


/iot.3«(/<)-i- »(/<,%  —  1  -F5) 
KW 


(1) 

(2) 


5  =  C(x30)  /C(x3h) is  the  coefficient  of  convective  diffusion. 

Hence,  the  problem  of  the  distribution  and  precipitation  of  a 
heavy  additive  (for  example,  in  the  presence  of  a  linear  constant  source), 
considering  both  the  gravitational  and  inertial  precipitation,  leads  to 
the  solution  of  the  diffusion  equation 


u  ( z ) 


<>C(x.z)  _  dC(x,  2) 

- - - -JJ  —  1  = 


dx 


dz 


-  2  - 


(3) 


with  boundary  condition  (1)  on  the  subjacent  surface1,  the  value  a, 
determined  by  Equation  2,  the  usual  condition  of  the  source  C(03z)  =■• 

Q 

=  (a  is  the  symbol  of  the  delta  function,  H  is  the  height 

u[tu 

of  the  source,  G  is  its  output),  the  usual  condition  at  infinity 
C  -*■  0  at  Jx2  +  x2  -*•  »,  and  isothermy  fAfaj  =  &a,  u(z)  =  u^a^). 

By  applying  the  Laplace  transform  for  a,  we  will  have 

*  '§  +  (*+ x)  4-^2— -i-  ’**»  <—«>.  (4) 


where  p  is  the  transformation  parameter  and  C(ptz)  is  the  transform  of 
the  function  C(x3z). 

The  solution  of  Equation  4,  obtained  by  the  method  of  variation  of 
the  derivatives  of  the  constants  and  correct  for  the  region  of  interest 
to  us  (a  <  //),  is 

C  (p.  z)  -  i4/C v  (a,  / p)  /,  (vx  Vp)  + 

At. ,  ifiVpw* /*>-/;(«* /?>) 

•  -T  M.  fc  Vri  lTf;^,7F?Twc^75)i (5) 

where  >1  =*  _=S_  f  £f/2t,  y,  =,  _L_  l/X  ,<»*)* , 

wnere  2  )  ’  1  H-?  K  *_ 

-,-£?/?  • 

F  „  /  ujklp-M  » 

—  ‘  a  • 

/  and  A  are  symbols  of  the  Bessel  functions  from  the  imaginary  argument 
and  the  MacDonald  function. 

The  solution  for  the  original  function  with  z  -  h,  obtained  by  the 
reverse  Laplace  transform,  is 

C  (a-,  h)  ---  C0  {x,  It)  Ci  (a’,  h), 


where 


Co  (A-,  A) - lf+q)kx  eXf)  L  ^  (H-  ?)s  A  j  /vL  *  (i -:-?)»*  J 


1  Condition  (1)  in  the  general  form  was  proposed  by  A.  S.  Monin;  we 
gave  it  a  concrete  form  when  applying  it  to  the  inertial  mechanism  of 
precipitation. 

-  3  - 


is  Rounds'  solution  [1]  for  2  =  h  for  the  same  problem,  but.  without- 
consideration  of  inertial  precipitation  (instead  of  condition  (1), 

condition  K(z )  — ^ — -*•  C  was  used,  with  2  -*■  0)); 

"  (x,  h)  = 

=  f  v Pi  (*  /p)  [g  /p  /w  (g»  Vp)  -  /,  (0*  Y~P)\  m 

"  0  Jco  ^VpW^VpJ  +  MajVp)  P‘ 


After  choosing  the  integration  curve  according  to  the  method  of 
M.  E.  Berliand  [2,  p.  27],  integral  (8)  can  be  presented  in  a  form 
suitable  for  numerical  integration  over  an  actual  variable  w: 


C  (X  m  \  2A  {Af  (FH  -f  QK)  -  /-  (Fl\  -  QH))  f**  , 

C x  (X,  It)  -  ^  A  +  (9) 


where 


r  «  ~  (AlA,  -  BxB*);  Q  «  ~  (AxB,  +  4A); 


H  —  E(oDi  -j-Cji  —  £aCx 

L  =  -f(£(oC3  -  A,);  M  =  ■—  (£coA  3  4-  B2); 

A  JJ i'?u)  — 7_v(vi(0)  _  D  _ 

~  sinitv/2  *  1  ~  cosstv/2  ’ 

„  i»(:>sM)— 7_v(»j(o)  e  Jv(t>3a)  fJ.,  foa)  # 

y]„  Sin  Vtv/2  *  2~  cosav/2 

„  .  D  _  +  ■ 

•™3  “  sinnv/2  ’  03  c6s"jtv/2 

Cx=  cos  (v  -{-  1 ) j  Jv+1  D\  —  sin  (v-f  1)]/ m  (fsa); 

C2=  cos  (21  j  7,  (vsa) ;  Z?2=  sin(21jjv  (t/sei); 

J  is  the  symbol  of  the  Bessel  function. 

The  density  of  the  precipitation  of  the  additive  on  the  ground  be¬ 
neath  the  vegetation  is 


go  =  %wc  {x,  h)\  (10) 

while  the  density  of  the  precipitation  of  the  additive  on  the  plants 
(relative  to  a  unit  area  of  the  ground)  is 
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(11) 


gp  =  \kh a  -f  w  (1  ~  £)1  C  (x,  h). 


For  practical  calculations  using  these  formulas,  it  is  necessary 
to  know  how  to  determine  the  new  dimensionless  characteristic  of  the 
cover  of  vegetation,  the  effective  entrapment  coefficient  a &h,  which 

d  wd2 

is  a  function  of  the  Stokes  criterion  Stk  =  m  ,  where  p  is  the 

18yU* 

particle  density,  Vy  is  the  viscosity  of  the  air,  and  X  is  the 
characteristic  size  of  the  barrier. 


The  direct  experimental  determination  of  this  characteristic, 
for  example,  by  the  direct  measurement  of  deposits  of  the  additive  on 
leaves,  is  laborious  and  difficult.  We  employed  an  indirect  means  of 
determining  a $h  =  /(Stk)  from  the  material,  balance  of  the  additive 
G  =  Gq  +  &  (Gq  =  total  amount  of  additive  deposited  on  the  ground  be¬ 


neath  the  plants,  G 


=  ditto,  on  the  plants): 


The  values  of  a &h  (taken  as  constants  in  the  range  Stk  >  8) ,  ob¬ 
tained  by  this  method  on  the  basis  of  the  results  of  field  tests,  as 
well  as  the  corresponding  values  ts3h}X  and  the  roughness  parameter  Zq 
are  listed  in  Table  2. 

Random  comparisons  of  the  experimental  data  with  solution  (6) 
were  conducted  with  these  values  of  the  parameters.  Agreement  was  found 
to  be  much  better  than  with  calculations  using  the  Rounds  formula.  Good 
agreement  between  the  calculated  and  experimental  values  for  tne  density 
gQ  on  the  ground  provides  a  basis  for  considering  that  the  density  of 

the  precipitation  on  the  vegetation  g  3  found  by  Formula  (11),  is  close 

to  the  actual  value. 


Calculations  show  that  consideration  of  the  inertial  precipitation 
by  the  root  method  changes  the  picture  of  the  process.  The  vegetation 
layer,  absorbing  the  additive,  causes  its  intense  diffusion  downward 
from  the  upper  layers,  i.e.,  as  if  it  were  sucking  the  additive  out  of 
the  atmosphere.  Consequently,  the  entire  field  of  concentrations  changes 
accordingly. 

The  exact  solution  (6)  given  above  has  a  value  in  principle,  but 
simpler  formulas  are  preferable  for  practical  approximations.  For  this 
purpose,  we  can  use  any  solution  of  the  diffusion  equation,  obtained 


5 


without  consideration  of  the  inertial  precipitation  (for  example. 
Rounds'  solution  Cq(XjZ)3  but  the  additive  current  must  be  divided 
into  two  parts  of  interest  to  us,  using  Equations  (10)  and  (lx) ;  from 
which  it  follows  that 


_  wC9  (x,0) 

“  1  -f  ti Zlagltti  (/i) ' 

(13) 

„  .  ^tiCo  (x,  0) 

(14) 

ap/m  (/i)/w  -}-  5  ‘ 

As  an  example,  in  Figure  1  we  have  plotted  the  results  of  deter¬ 
mining  gQ  according  to  the  exact  formulas  (6),  (7)  and  (9)  which  we 

have  found  (Line  1),  according  to  Rounds'  formula  (7)  (2),  and  accor¬ 
ding  to  approximate  formula  (14)  with  determination  of  Cq(,Xj  0)  by 

Rounds'  formula  (7)  (3).  In  the  calculations,  we  used  the  following 
values  of  the  parameters:  G  =  802  mg/m,  tf  =  1.6  m,  h  =  0.2  m, 
k  =  0.166  m/sec,  w  =  0.0208  m/sec,  q  =  0.281,  u ^  =  4.12  m°*71Vsec, 

5  =  1.0,  a &h  =  0.486,  u(h)  =  2.62  m/sec,  K(h)  =  0.0332  m2/sec.  The  cor¬ 
responding  experimental  values  of  gQ  are  represented  by  points. 


Table  2 


1 

J  *8,  CM 

1 

5 

ft,  CM 

ajft 

X,  CM 

Steppe,  natural 

0,52 

1,0 

5 

0,0015 

— 

vegetation  cover 

0,86 

1,8 

1,0 

1,0 

8 

20 

0,006 

0,14 

0,20 

2,7 

1,0 

30 

0,055 

0,3S 

Wheat  in  blossom, 
50-80  plants/m2 

0,38 

60 

0,20 

0,30 

It  is  clear  from  the  graph  that  both  the  exact  solution  and  the 
approximate  formula  (14)  are  much  closer  to  the  experimental  data  than 
Rounus'  solution;  as  we  have  already  pointed  out,  this  provides  a  basis 
for  considering  that  formula  (13)  is  reliable;  a  direct  experimental 
verification  of  this  formula  would  be  laborious. 

The  possibility  of  calculating  separately  the  precipitation  current 
on  the  ground  below  plants  (loss  of  chemicals)  and  the  precipitation 
current  on  the  plants  (amount  of  chemicals  employed  to  advantage)  is  es¬ 
pecially  important  for  agricultural  applications. 
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